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1. Introduction 

Explicit sharp estimates for the Green function of the Laplacian in C^'^ domains 
were completed in 1986 by Zhao [43]. Sharp estimates of the Green function of 
Lipschitz domains were given in 2000 by Bogdan [6] . Explicit qualitatively sharp 
estimates for the classical heat kernel in C^'^ domains were established in 2002 
by Zhang [42]. Qualitatively sharp heat kernel estimates in Lipschitz domains 
were given in 2003 by Varopulous [41]. The development of the boundary po- 
tential theory of the fractional Laplacian follows a parallel path. Green function 
estimates were obtained in 1997 and 1998 by Kulczycki [29] and Chen and Song 
[21] for C^'^ domains, and in 2002 by Jakubowski for Lipschitz domains [28]. In 
2008 Chen, Kim and Song [19] gave sharp explicit estimates for the heat kernel 
Poit, X, y) of the fractional Laplacian on C^'^ domains D. The main contribution 
of the present paper is the following result. 

Theorem 1. If D is n-fat then there is C = C{a,D) such that 

c~^p^{td > t)py{TD >t)^ PDij^^^y) ^ cp^TD > t)py{TD > t) (i) 

for < t ^ 1 and x,y ^ D. 
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Here x, y) is the heat kernel of the fractional Laplacian on M"*, and 

P^(ri) >t)= pD{t,x,y)dy 

defines the survival probability of the corresponding isotropic a-stable Levy pro- 
cess in D. The result applies also to unbounded domains, in particular to do- 
mains above the graph of a Lipschitz function, where we can take arbitrary 
t > 0. In fact (1) holds with C = C(a, d, k) under the mere condition that D is 
(K,t^/")-fat at X and at y, see Section 3 and 4 for definitions and results. For 
exterior domains wc have a result free from local geometric assumptions: 

Corollary 1. //diam(£''^) < oo then (1) holds with C — C{a,d) for all t > 
diam(D'^)" and x,y D. 

For exterior domains of class C^'^ a more explicit estimate is given in Theo- 
rem 3 below. We also like to note that a useful variant of Theorem 1 is given in 
Theorem 2. 

(1) is motivated by these applications of the semigroup property oi po- 



PD{2t,x,y)= / pD{t,x,z)pD{t,z,y)dz ii, P'^iTD > t)c{t) , 
where c{t) = sup^ y^^dp{t,z,y) ^ sup^ y^g^a poit, z,y), see (12), and 



PD{3t,x,y) = J J pDit,x, z)pd it, z, w)pd {t,w, y)dwdz 

^ P''{TD>t)c{t)py{TD>t). 

The latter inequality is quite satisfactory for x — y, because c{t) = p{t,x,x). 
Off-diagonal (x, y) in (1) require, however, a deeper analysis. Our proof of (1) 
is based on the boundary Harnack principle (BHP) [14] (see also earlier [40]), a 
version of the Ikeda-Watanabe formula (18) ([27]), scaling (14), and comparabil- 
ity of p with its Levy measure (5), see (28). Counterparts of these are important 
in view of possible generalizations. 

In what follows (1) and analogous sharp estimates will be written as 

PD{t,x,y) ^ P'{td > t)p{t,x,y)py{TD > t) , 

meaning that either ratio of the sides is bounded by a number C € (0, oo), and 
C does not depend on the variables shown (here: t, x, y). We will skip C from 
notation if unimportant for our goals. 

Let 5d{x) = dist(a;, D'^). As mentioned above, domains D of class C^^^ enjoy 
the following sharp and explicit estimate of Chen, Kim and Song [19]: 




PD{t,x,y ) 
Pit,x,y) 



<t i^l, x,y eM.'^ . (2) 
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Wc note that (2) agrees with (1) because by [10, CoroUary 1], 

In fact, starting with (1), we are able to recover and strengthen (2), with a 
simpler proof, see Example 5 and Proposition 1 below. We note that (1) was 
conjectured in [10] based on the cases of C^'^ domains [19] and circular cones 
[10]. We should also mention that the Gaussian estimates of Varopoulous [41] 
have a shape similar to (1), in particular they involve the survival probability. 
Thus the present paper builds on the evidence accumulated in [41], [19] and [10]. 
Wc also note that the upper bound in (2) was proved in 2006 by Siudeja for 
semibounded convex domains [39, Theorem 1.6] and stated for general convex 
domain in [39, Remark 1.7]. Some of our present techniques were inspired by 
[32, Theorem 4.2] of Kulczycki and Siudeja, [2. Proposition 2.9] of Bafiuelos and 
Kulczycki, and [1, Section 4] of Bogdan and Bafiuelos. 

It is a consequence of Lemma 1 below that we can apply BHP [40, 14] to 
conveniently estimate P^{to > 1) by some kernel functions of D, namely by 
the Martin kernel with the pole at infinity or the expected survival time (we 
use scaling to estimate P^{td > t) for general t > 0). The estimate and the 
resulting bounds for the heat kernel are collected in Theorem 2, followed by a 
number of applications. In particular we give a simple proof of the main result 
of [10] for the circular cones V: 



pv{t,x,y) (.I^ITT^J (.I^TTatJ 



(3) 



Here f3 £ [0, a) is a characteristics of the cone, and all t > and x,y gM.'^ are 
allowed. We should add to (1), (2) and (3) that ([4], [16]) 

t>0,xem'. (4) 

Here c = c(a, d), meaning that c e (0, oo) may be so chosen to depend only on 
d and a. We like to note that the estimates for general K-fat domains cannot 
be as explicit as those for C^'^ domains. In particular the decay rate (3 at 
the vertex of a cone delicately depends on the aperture of the cone, sec [10, 
1, 35] (see also [6]). Nevertheless Lipschitz domains offer a natural setting for 
studying the boundary behavior of the Green function and the heat kernel for 
both the Brownian motion and the isotropic a-stable Levy processes. This is 
due to the scaling, the rich range of asymptotic behaviors depending on the 
local geometry of the domain's boundary, connections to the boundary Harnack 
principle and approximate factorization of the Green function, and applications 
in the perturbation theory of generators, in particular via the 3G Theorem [6, 
1; 43, 26, 7] and 3P Theorem [13]. The K-fat sets are a convenient generalization 
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of Lipschitz domains, with similar features. Noteworthy, (1) is an approximate 
factorization of the heat kernel (see also [6, 14] in this connection). 

We should add that the C^'^ condition specifies the geometry of a domain 
only in bounded scales (see Definition 3). This renders the range of time in (2) 
restricted to < t ^ 1. In what follows we will also study the probability of 
survival for large times (and unbounded domains). This is straightforward for 
special Lipschitz domains (thus for circular cones), but less so for general K-fat or 
C^'^ domains. As an interesting case study we consider domains with bounded 
complement (i.e. exterior domains) of class C^'^. These have distinctive geome- 
tries at infinity and at the boundary, resulting in nontrivial completion of (2). 
We remark that exterior C^'^ domains in dimension d > a have been recently 
studied in [22], too. We also remark that [25, Theorem 4.4] bounds the survival 
probability of the relativistic process in a half-line, and [31] gives an explicit 
formula for the transition density of the killed Cauchy process (a = 1) on the 
half-line. Regarding other recent estimates [20, 17, 36, 23, 3] for the transition 
density and potential kernel of jump-type processes, we need to point out that 
generally these only concern processes without killing. Killing corresponds to 
the Dirichlet "boundary" condition (analogous to negative Schrodinger pertur- 
bation [8, 12]) and it severely influences the asymptotics of the transition density 
and Green function. Needless to say, the asymptotics is crucial for solving the 
Dirichlet problem [24, 25]. 

We like to mention possible applications and further directions of research. 
The estimate (1) fits well into the technique of Schrodinger perturbations of 
[12], which should produce straightforward consequences. Also, the distribution 
of td, given by (18) below, can be estimated by using (1). Further, we conjecture 
that for certain domains D, \iva.pD{t,x,y)/ P^(td > t) exists as x approaches 
a boundary point of D. This may lead to representation theorems for nonneg- 
ative parabolic functions of the fractional Laplacian (compare [14, Theorem 2 
and 3]), and construction of excursion laws. We need to remark here that our 
estimates are inconclusive about the (irregular [14]) boundary points of D, but 
we conjecture that (1) indeed extends to dD. Finally, it seems important to 
understand the behavior of pD{t,x,y) for domains which are rather small at a 
boundary point or at infinity. In this connection we refer the interested reader 
to the recent study of intrinsic ultracontractivity by Kwasnicki [33], see also 
[30, 19, 10] and the notion of inaccessibility in [14]. 

Our general references to the boundary potential theory of the fractional 
Laplacian are [7] and [14]. We also refer the reader to [9] for a broad non- 
technical overview of the methods and goals of the theory. 

The paper is composed as follows. In Section 2 we recall basic facts about the 
killed isotropic a-stable Levy processes. In Section 3 we prove Theorem 1 and 
Corollary 1. In Section 4 we state and prove Theorem 2 and give applications to 
specific domains. In particular we strengthen (2) and part of the results of [19] 
(see Proposition 1, Theorem 3 and Corollary 2), and we discuss exterior C^'^ 
domains in dimension d = I < a. 
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2. Preliminaries 

In what follows, denotes the Euclidean space of dimension d ^ 1, is the 
Lebesgue measure on M**, and < a < 2. Our primary analytic data are: a 
nonempty open set D C R** and the Levy measure given by density function 

2"r((rf+a)/2) ,^,_,,, 

Ky) - ,./2|r(_«/2)| 1^1 ■ 

The coefficient in (5) is such that 

/ [1 - cos(e • y)] i.(y)d2/ = , C G R'^ . (6) 
For (smooth compactly supported) (p £ C^(R'^), the fractional Laplacian is 



A"/20(2;) ^ lim / [0(a; + y) - (l){x)] v{y)dy , x- G R'' 

eiO J|y|>e 

(see [7, 9] for a broader setup). If r > and (f)r{x) = 4>{rx) then 

A"/2(/.^(a;) =r"A"/2(/.(ra;), x G M'^ . (7) 
We let pt be the smooth real-valued function on R'' with Fourier transform, 

/ ft(a;)e" «dx = e-*l«l°, t > , ^ e R"* ■ (8) 

In particular the maximum of pt is pt(0) = 2^~°'TT-'^/'^a-^T{d/ a) /T{d/2)t~'^/°' . 
According to (6) and the Levy-Khinchine formula, {pt} is a probabilistic convo- 
lution semigroup with Levy measure v{y)dy, see [38], [16], or [9]. We have the 
following scaling property, 

p^{x) = <-''/>i(i-i/«x) , t>{),xeW^, (9) 

which may be considered a consequence of (8). Noteworthy, by (4) we have 

Pt{x) ^ P2t{x) , i>0,a:GR''. (10) 

We denote 

P{t,x,y) =pt{y-x), 

and we have 

oo 

J jp{u — s,x,z) 9ti(/)(u, z) -t- A"^^(^(u, z) dzdu — —4>{s,x) , (11) 
where s G R, x G R'', and cf) G C~(R x R''), see, e.g., [12, (36)] . 
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We define the isotropic a-stable Levy process {Xt,P^) by stipulating transi- 
tion probabifity 

Ptix,A) = / p{t,x,y)dy, t > , x eR'^ , A cR"^ , 



J A 

initial distribution P^{X{0) = x) = 1, and cadlag paths. Thus, , denote 
the distribution and expectation for the process starting at x. We define the 
time of the first exit from or survival time, 

td = inf{t >{): XtiD}, 

and the time of first hitting D, 

Td = mi{t >0: XteD}. 

We define, as usual, 

PDit,x,y) =p{t,x,y) - E^'Itd < t; p{t - TD,Xr^,y)], t > 0, x,y e M.'^ . 
We have that 

< pnit, X, y) = pD{t, y, x) < p{t, x, y) , (12) 

hence 

j PD{t, X, y)dy = j poit, x, y)dx ^ 1 . (13) 

If a; G D'^ is regular for the Dirichlet problem on D [14], i.e. P^{to = 0) = 1, then 
PD{t, x,y) = and (1) is trivially satisfied. By this remark, if all the points oidD 
are regular for D, then we can write x,y G R'' in Theorem 1, instead of a;, y £ D. 
The remark also applies to Examples 1-8 in Section 4. By the strong Markov 
property, is the transition density of the isotropic stable process killed when 
leaving D, meaning that we have the following Chapman-Kolmogorov equation, 

Pd{s,x, z)pD{t, z,y)dz = pd{s + t,x,y) , s, t > 0, x, y S M'' , 

and for nonnegative or bounded (Borel) functions / : M'' ^ M, 

f{y)pDit,x,y)dy^E'[TD<t; /(X*)] , t > , x e . 

For s e R, a; g R'', and (j) € Cf'(M x D), we have 

Pd{u — s,x,z) 9„0(it, z) -f A"^^0(m, z) dzdu = —(t){s, x) , 



s D 



which extends (11) and justifies calling pn the heat kernel of the (Dirich- 
let) fractional Laplacian on D. It is well-known that is jointly continu- 
ous and positive for {t,x,y) £ (0,oo) x D x D. We have a scaling property: 
Proir^t, rx, ry) = r'^pnit, x, y), r > 0, or 

PD{t, X, y) = t-'*/>,_i/„^(l, t-i/"y) , i > , x, J/ e R^ (14) 
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in agreement with (9) and (7). Thus, P^'^iTro > = ^"^(td > t), or 

P^rD>t)= [ pD{t,x,y)dy = P''""'-{T,-u.D>^)- (15) 

Remark 1. For c > consider v ~ cv, the corresponding heat kernels p, po, 
probability and expectation P^, . Clearly, pu{t,x,y) = pu{ct,x,y). 

The Green function of D is defined as 

GDix,y)= / pDit,x,y)dt, (16) 
Jo 

and scaling of p^i yields the following scaling of G^j , 

GrD{rx,ry)=r"-''GD{x,y). (17) 

A result of Ikeda and Watanabe [27] asserts that for x S 13 the P^-distribution 
of {TD,Xro^,Xro) restricted to Xr^-^ ^ Xt-^ is given by the density function 

(s, It, z) H> p£)(s, X, m)i'(z — u) . (18) 

For geometrically nice domains, e.g. for the ball, P^(Xt-q_ ^ X^jj) = for 
x G D [14], and then by (16) and (18) the P^-distribution of Xr^ has the 
density function given by the Poisson kernel: 

Pd{x,z)^ / Gd{x,u)i'{z — u)du . (19) 
Jd 

For Xq G M'^ and r > we consider the ball B{xo,r) = {.t £ R'^ : [a; — a;o| < r} 
and B'^{xo,r) = {x : \x — xo\ > r} (open complement of a ball). 
There is a constant C depending only on d, a and p, such that 

(J 

Pu{xi,yi)Pu{x2,y2) ~ Pu{xi,y2)Pu{x2,yi) , (20) 

whenever U C B{xo,r) C R'^ is open, < p < I, r > 0, xq € R'^, xi,X2 S 
U n B{xo,rp), and yi,j/2 G B(xQ,rY. This boundary Harnack principle (BHP) 
follows from [14, Lemma 7 and the proof of Theorem 1], and it is essentially an 
approximate factorization of Pjj- We encourage the interested reader to directly 
verify the estimate in the special case of (22) below. 

The Green function and Poisson kernel of B{xq, r) are known explicitly: 

a/2-1 

GBi.„r)ix,v)^Ba,^\x-vr'' (21) 



p f \ r \r^-\x~XQ\^ 

PB{xo,r)(X,y) = Cd 



y — xo\'^ — 



a/2 ^ 

' ^^^^ 



where Bd,a = r(d/2)/(2"7r''/2[r(a/2)]2), C^,, = r(d/2)7r-i-'^/2 sin(™/2), 

^ = (^2 _ 1^ _ ^^|2)(^2 _ _ ^^|2)/|^ _ ^|2 ^ 
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\x — Xq\ < r, \v ~ xo\ < r, and \y — xq\ ^ r; sec [5], [37]. Thus, 
P^\Xr,^„J> R) = / PB(o,i)(.x,2/)dy « ^ 

J\y\:^R ^ 

where x € -8(0, 1) and 2. Also, for |a; — xo\ ^ r we have ([8]) 



(23) 



All the sets and functions considered below are Borelian. Positive means 
strictly positive. Domain means a nonempty open set (connectedness need not 
be assumed in this theory). 



3. Factorization 

We consider nonempty open set _D C M'^. 

Definition 1. Let x ^ D, r > and < k ^ 1. We say that D is (K,r)-fat at 
X if there is a ball B{A, kt) C D D B{x, r). If this is true for every x G D then 
we say that D is (k, r)-fat. We say that D is K-fat if there is i? > such that D 

is (k, r)-fat for all r G (0, R]. 

Remark 2. The ball is 1/2-fat. 




Definition 2. Given B{A, k) C D n B{x, 1), we consider U = D n B{x, \x - 
A\ + k/3), Bi ^ B{A, k/'S) c U and B2 = B{A\ k/6) such that B{A', k/'S) C 
B{A, k) \ U, see the picture. 

Lemma 1. There is C = C'{a,d,K) such that if D is (k, l)-fat at x then 

P"(rz5 > 1/3) 5:CP^(Tp >3). (25) 
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Proof. Consider x E D and B{A, k) and U as above. For |a; — A| < k/2, 

> P^{rB(..,./2) > 3) = P°(rB(o,,/2) > 3) > , 

and (25) is proved. We will now assume that |x — A| ^ k/2. Wc note that 

P^TD > 1/3) s; P'^Tu > 1/3) + P'^Xru e D) . (26) 

We have P^^iX^^, e D) ^ Pu{x,y)dy. Indeed, li B = B{x, \x - A\ + k/3) 
as in the definition of f/, then P^(X^„ € dU n D) ^ P^iXr^ e dB) = 0, see 
the discussion preceeding (19) above. Similarly, P^{Xr^, G B2) is an integral of 
the Poisson kernel Pjj- We consider BHP for xi = x, X2 = A, p = 1 — k/3 > 
(1 — k)/(1 — K + k/3). Integrating (20) on D and B2 we obtain 

P^iXru e ^) " ""p^iXru e B2) ■ 

Wc note that (the denominator) P^^iXr^ e B2) ^ P'^iXm, G -B2) ^ c > 0, 
see (22), therefore P'=(X^„ G £») < cP"=(X^„ G 52)- We also observe that 
u i-> Jg^ v[y — u)dy is bounded away from zero and infinity on U . By (19), 



P^{Xr„<^B2)= i Gu{x,u) I v{y — u)dyduKi / Guix,u)du = E^tjj . 
Ju J B2 Ju 

Clearly, P^(t[7 > 1/3) ZE^'tu. By (26), P^(ti5 > 1/3) cE='tu. By the 
strong Markov property, 

E-Tu s; cP^{Xr„ G P2) ^ cE- [Xru e P2; P""-- (tb(X.^ ,k/6) > 3)" 
^ cP^(rz,>3). 

□ 

Remark 3. If I? is (k, l)-fat at .x then by the above proof wc have 

P^TD > 1/3) « P^TD > 3) « P"(tz, > 1) « P^Xr, eD)^ E-Tu . (27) 

In fact, we can replace 3 by any finite £ ^ 1, at the expense of having the 
comparability between each pair of expressions in (27) holding with a constant 
C = C{a, d, K, £). 

Lemma 2. Consider open Di,D^ C D such that dist(_Di, Da) > 0. Let D2 = 
D\{DiU D3). Ifx G Di and y G D3, then 

PD{^,x,y) i^P'^iXro^e D2) sup p{s,z,y) + E'^td^ sup v{z - u) , 

s<i,zeD2 ueDi.zeDs 



and 



PD{l,x,y)^P^TD,>l)py{TD,>l) inf iy{z~u). 

ueDi , zeDs 
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Proof. By the strong Markov property 

PD{l,X,y) = S^bnll - TDi,Xro^,y),TDi < 1], 

which is 

E''[pDil-TD„Xr^^,y),TD,<l,Xr^^eD3] = I + II . 

Clearly, 

l^P'^iXr^^ eD2) sup pis,z,y). 

S<1, z£D2 

Consider Di such that P^{Xt^ £ dDi n Z?) = 0, e.g. Di being an intersection 
of D with a Lipschitz domain. By (18), the density function of (r^i^ , X^-^^^ ) at 
(s, z) for z £ D equals 



f^{s,z)= j pDi{s,x,u)i'{z — u)du . 

For z e D3, 



f^{s,z)^ / pjj^{s,x,u)v{z — u)du ^ P^ {t£,^ > s) sup v{z — u) , 

J Di ueDi, zeDs 

hence, by (13), 
II = [ [ poil- s,z,y)f''{s,z)dzds 

JQ J D3 

^ sup i'{z ~ u) / / pd{1 — s, z,y)P^ {td-^ > s)dzds 

ueDi,z£D3 Jq Jds 

^ / P^{tdi > s)ds sup v{z ~ u) ^ E''' Tjj^ sup v{z — u) . 
Jo u<£Di,zeDi ueDi,zeD3 

The upper bound follows. The case of general Di follows by approximating from 
below, and continuity of p and v. The lower bound obtains analogously: 

II ^ inf i/(z - m) / / poil - s,z,y)P'^{TDi > s)dzds 
ueDuzeDs Jo Jd^ 

> P^{td, > 1) inf ,^{z-u) r I pD3{l 
ueDi.zeDs Jq Jij 



s, z, y)dzds . 



□ 



Remark 4. Lemma 2 also holds for i), p, P^ and i?^ of Remark 1. 
In what follows we will often use the fact that 

lAi^(z-u) «p(l,u,z). (28) 
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Lemma 3. If D is (k, l)-fat at x and y then 

Pd{2, X, y) ^ Cia, d, P'' {td > 2)P^(rz5 > 2)p(2, y) . 

Proof. If \x — y\ ^8 then p{l,x,y) « 1, and by the semigroup property, (10) 
and Lemma 1, 



PD{l,x,y) = I pD{l/2,x,z)pD{l/2,z,y)dz 
^ supp{l/2,z,y)P%TD>l/2) 

Z 

^ cP^{TD>l)p{l,x,y). (29) 

Here c = c(a, d, k). If |x — ?/| > 8 then we will apply Lemma 2 with Di = 
U = Dn B{A, \x- A\+ k/3), as in Definition 2, and D3 = {z e D : \z - x\ > 
\x-y\/2}. Since sup^<i_ ^^^^ p{s, z, y) ^ cp{l,x, y), and sup^^^,^^ ^g^,^ i^{z-u) < 
cp(l,x,y), see (28), by Remark 3 we obtain, 

PD{l,x,y) < cpil,x,y) [P"(X,„ £ D) + E^tu] ^ cP^td > 1) p(l, x, y) , (30) 

hence by (29), (30), symmetry, semigroup property and Lemma 1, 



Pd(2,x,2/) = / pDil,x,z)pDil,z,y)dz 



^ cP%TD > l)py{TD >l) J p{l,x,z)pil,z,y)dz 
< cP^TD > 2)py{TD > 2)p{2, X, y) . 



□ 



Under the assumptions of Lemma 3, C = C{a, d, n) exists such that 

PD{l.x,y)^CP^TD > 1)P"(tz5 > l)p(l,x,y). (31) 

Indeed, according to Remark 1 we consider v — and the corresponding p, 
p_D, P^, obtaining 

PD{l,x,y) = pD{2,x,y)^CP^[TD>2)P^TD>2)p{2,x,y) 
= CP^{TD>l)P^{TD>l)p{l,x,y). 

Lemma 4. If r > then there is a constant C ~ C{a,d,r) such that 

PB{u,r)UB{v,r){i,U,v) ^ Cp{l,U,v), U,veM.'^. 

Proof. For |u — w| ^ r/2 we use (28) and Lemma 2 with D = B{u, r) U B{v, r), 
Di = B{u, r/8) and = B{v, r/8): 

PB(u.r)uB(v.r){i,u,v) ^ P"(rD^ > 1)P''(ti33 > 1) iuf iy{z - u) 



c[^'°(TB(0,r/8) > 1)]%(1,M,V). 
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For |u — w| ^ 7'/2 by (4) wc simply have, 

PS(«,r)UB(t),r)(l,'", W) > , inf , PB(o,r) (1 , 0, z) ^ C ^ Cp(l , W, w) . 

\z\<r/2 

□ 

Lemma 5. If D is (k, l)-fat at x and y then 

PDi3, X, y) ^ C{a, d, ti) P^'iro > 3)P^(rc > 3)p(3, x, y) . 

Proof. Consider t/^, , and U^, selected according to Definition 2 for x 
and y, correspondingly. By the semigroup property, Lemma 4 with r = k/6, and 
(4), 

PD{3,x,y) ^ pD{l,x,u)pD{l,u,v)pD{l,v,y)dudv 
J 3% Jb? 



cp{\,x,y) I pD{l,x,u)du / pD{l,v,y)di 
Jbs JbH 



For u(^B^= B{A\ k/6), by Lemma 2 with Di^U"" = U and A3 = B{A' , k/4), 
and by Remark 3 we obtain 

Pd{1,x,u) > P'-'iru > 1)P^{tb(o.k/v2)> I) ini v{z - w) 

weu.zeDs 

^ cP^Tu >1)^cP''{td>1). 
Similarly, piy{l,v,y) ^ cPy{TD > 1), hence, by Lemma 1, we have 

PD{5,x,y) ^ cP\Tu>\)p{:i.x,y)P%Tu>\) 
^ cP\td > i)p{i,x,y)P^{TD > 3) . 

□ 

Under the assumptions of Lemma 5 we also have that 

PDil,x,y)^Cia,d,t,)P''iTD > 1)P^(td > l)pil,x,y). (32) 
This is proved analogously to (31). 

Proof of Theorem 1. Assume that R ^ 1 and D is (K,r)-fat for < r ^ i?. If 
t^/" e {0,R] then <-i/"£) is {k, l)-fat. The estimate (1) follows from (31), (32), 
and scaling, see (14) and (15). In fact we have C = C{a, d, k) in (1). If i? < 1 
then we argue as in the case of (31) C = C{a, d, k, R), or, alternatively, we use 
Remark 6 below. □ 

Proof of Corollary 1. Note that D is (l/4,r)-fat for r ^ 2diam(Z?'^). and so 
we obtain (1) for t ^ 2"diam(D'^) with the same constant C. If we consider 
v = 2~°'v and argue like in the case of (31), then we obtain the wider range of 
as in the statement of Corollary 1. □ 
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Remark 5. Since the k- fatness condition is more restrictive when k is bigger, 
the above constants C = C{a,d,n) may be chosen decreasing with respect to 
K. Also, if D has a tangent inner ball of radius 1 at every boundary point then 
the constants in Lemma 3 and 5 depend only on a and d. 

Remark 6. If D is (k, r)-fat at x and 1 ^ K < oo then D is [n/K, rK)-ia.i at x. 
This observation together with scaling allows to easily increase time, compare 
(31) or (32), at the expense of enlarging the constants of comparability. The 
argument, however, does not allow to decrease time. Remark 1 is more flexible 
in this respect. 



4. Applications 

Wc let sd{x) = Ej-td = J GD{x,y)dv if this expectation is finite for x € D, 
otherwise we let s jj (x) = Mjj (x) , the Martin kernel with the pole at infinity for 

sr t \ V Gpjx^y) 
Md{x) = lim — — . 

We should note that this (alternative) definition of S£) is natural in view of [14, 
Theorem 2]. The choice xq G D is merely a normalization, Md{xo) — 1, and 
will not be refiected in the notation. By the scaling of the Green function (17) 
we obtain 

Sroirx) sd{x) 

r^ = — K' x,yeD,r>0. 33) 

SrD(ry) SD{yj 

We denote by Ar{x) or Ar{x,K,D) every point A such that B{A,Kr) C D Ci 
B(x,r), as in Definition 1. Noteworthy, Ar{x) approximately dominates x in 
terms of the distance to dD: 

SD{Ar{x))^rVSD{x). (34) 

If D is (k, l)-fat at x then rD is (/t, r)-fat at rx, and (every) rAi(x, n, D) may 
serve as Ar{rx, n, rD). 

Theorem 2. If D is (K,t^/")-fat at X and y then 

P''{rD>t)^ .. , (35) 

SD[Ati/^{x)) 

where C = C{d,a,K) and, furthermore, 

PD{t,x,y)K — 1-. i-^p{t,x,y) — — (36) 

SD{Ati/^{x)) SD(Ati/„(y)) 

Proof. To verify (35) we first let i = 1 and assume that D is (k, l)-fat at x. Let 
A = Ai{x). If E^To < oo then we consider the set U C D oi Definition 2, and 
we obtain 

E^'td^ E'^Tu + E^'soiXr^). 
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By Remark 3, E^tu ~ P^{td > !)• Since E^tu w 1, we trivially have 
Similarly, P^{Xru G -D) w 1. By BHP and Remark 3 we obtain 

E'-'soiXr^) P"(Xr„ e D) 



P^'-iXr,, eD)^P-irD>l). (37) 



This yields (35) in the considered case. If E^to = oo then is harmonic and 
we have sd{x) = E^SDiXr^) (see [14, Theorem 2 and (77)]) and we proceed 
directly via (37). The case of general t in (35) is obtained by the scaling (33) and 
(15). Finally, (36) follows from (35) and Theorem 1. The resulting comparability 
constants depend only on a, d and k. 

□ 

Remark 7. Assume that D is K-fat, so that there is i? > such that D is (k, r)- 
fat for every r ^ R. Then (35) and (36) hold with C = C{d, a, k) for all x,y £ D 
and t ^ R". 

Below we give a number of applications. 

Example 1. We let i? > and D = B{0, R) C M'^. By (24), the expected sur- 
vival time is sd{x) w S'^^^{x)R"^^, where C = C(d,a). By (34), 5^5(^41/0 (x)) w 
(i^/" V 6Dix))"/'^R"/'^, therefore for ah t i?" and x,y G M'^, 

f-fa>Og .,/:?T' . -flA%fe^)°''. (38) 



(tV" v<Si5(x))"/2 V 
and 

PD{t,x,y)S,(^lA^-^y{t,x,y)(^lA^-^^ , (39) 

To be explicit, 6b{o,r){x) = {R— l^^l) VO, and 6b{o,r)''{x) = {\x\ — i?) VO, and 
(38), (39) on D'^ follow because all x € D'^ are regular for D. 

Example 2. Let C M'' be a halfspace. The Martin kernel with the pole at 
infinity for D is sd{x) = 5'^'^{x) ([1]). We see that (38) and (39) hold with 
C = C{d, a) for all t S (0, 00) and x,y€W^. 

Example 3. Let D = B'=(0, 1) C and d ^ a. By Kelvin transform ([18] or 
[14]) and (21), 

N"'1yr''GB(^,^) _ |x|°-'^Gb(^,0) 
^^""^ - A-|xo|"-'^|2/|«-^Gb(^,^) " |xo|"-'^Gs(^,0)' 

where 

/-kr^-l „a/2-l 
GB(z,0)=Brf.„|z|"-M ^^^ds, 0<|z|<l. 
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Thus, there is c = c(a;o, a) such that 

If d > a then sd{x) w 1 A5^''^(x), sc(A4i/„ (x)) w 1 A (t^/" V (5d(x-))"/2^ thus 
^ ^ ^ 1 A V 5d (x))"/2 (1 a <i/")«/2 ' 

and 

for all < t < oo and x, ?/ € M''. Here C = C{d, a). 

For a = d = 1, (40) yields sd{x) « log(l + 5]^\x)), SD{Ati/.{x)) « log(l + 
{t V (5£)(a;))^/2) thus for ah < t < oo and x, y e R'' wc have 

^ > ^ log(l + (i V <5z,(a:))i/2) - ^ log(l + tl/2) ' ^^^^ 

and 

PD{t,x,y) 
p{t,x,y) 

Sharp explicit estimates for Pb'-{o.r) with arbitrary R > follow by scaling. 
Example 4. Let D = B=(0, 1) C R'^ and 1 = d < a. We have that 

G{o}c(x,y) = Gz3(a;,y) +£;'^G{o}c(^TB,y) • 
Let Ca = [-2r(a) cos(7ra/2)]^^ By [18, Lemma 4] for x,y eR, 

G^or i^, y) = Ca (iyr-' + kr-i - |y - . 

If follows that 

Goix^y) = Co {\xr~' -\x~ yr-i - E%\Xr^r-' - \X,^ - y^-') . 

Since \Xr^\ ^ 1 a.s., \iuiy^oo{-\x - yf-^ + E''\Xr„ - = 0, for every 

a; £ R. If |a;| > 2 then we can find c = c(q;, a;o) such that 



1 A 



log(l + j£(x)) \ / log(l + ^V2(zy)) \ 

l0g(l+<l/2) J \^ l0g(l+il/2) ) 



md{x) = Ni° ^ ^""'^""i" ' =0(1x1"-^ ~^"|x.„r-i) 

|.Tr-i «<5d(.t)"-i 



On the other hand, by BHP Miy{x) « 5^'^{x) ii5(x) ^ 1 (compare Example 2). 

We thus have sd{x) « 5d~^{x) A 5'^^{x), sd(Ai/»(x)) « (t^/" V fo(x))"~i A 
(*!/" V i5D(a;))"/2, and for &\\ Q < t < oo, x,y € W\ we obtain 

P-(r , .^g. ^r^(x)A^g,/^x) 

^ (tV"vfe(x))«-iA(tV"vfo(x))"/2 ' ^ ^ 
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hence 

Here C ~ C{a). To estimate Pb':{o,b.) with arbitrary i? > we use scahng. 

Definition 3. We say that (open) D is of class C^'^ at scale r > if for every 
Q edD there exist balls B{x',r) C D and B{x" ,r) C Z?'' tangent at Q. If £> is 
C^'"^ at some (unspecified) positive scale (hence also at smaller scales) then we 
simply say D is C^'^. 

C^-^ domains may be equivalently defined using local coordinates ([34]). 

Remark 8. If D is C^'^ at scale r then it is (l/2,p)-fat for all p G (0,r]. 

Remark 9. Let be C^'^ at scale r. Let x ^ D, and let Q e 9Z3 be such that 
Sd{x) = \x — Q\. Consider the above balls B{x',r) and B{x",r). If Soix) < r 
then let B^ ~ B{x',r), otherwise B^. = B{x, 6d{x)). Thus Sb^{x) = 6d{x), and 
the radius of B^; is r V Sd{x). 

Example 5. We will verify (2) for C^'^ domains D. For the proof we initially 
assume that D 7^ R'' is C^'^ at scale r = 1. Let x € D. Wc adopt the notation 
of Remark 9 and consider (the ball) B^ and (the open complement of a ball) 
B=(x", 1) tangent a.t Q G dD. Since B^ C D C B=(.t", 1), wc have 

P^iTB^ > 1) ^ P^TD > 1) ^ P"(tbo(,„,i) > 1) . 

Clearly, SbA^) = Snix) ^ \Q - x\ = dB^(.",i)(2:)- By (38) and (41)-(43), 
By Remark 8 and Theorem 1 there is C = C{d, a) such that for all x,y G M'^, 

(,.M^)(_i^.,-..)(,.M|f!), 

If D is C^'^ at a scale r < 1 then r'^D is C^'^ at scale 1. This yields (2) in 
time range < t ^ r". Remark 3 allows for an extension to all t G (0, 1], with 
a constant depending on d, a and r. The case of D = R'' is trivial. 

Further estimates for C^'^ domains will be given in Proposition 1, Theorem 3 
and Corollary 2. 

Example 6. Let d > 2. For x = {xi, . . . , Xd-i,Xd) G R'' we denote x = 
{xi, . . . , Xd-i), so that X — (i, x^). Let A < 00. We consider a Lipschitz func- 
tion 7 : R''^^ R i.e. \^{x) — 7(y)| ^ X\x — y\. We define a special Lipschitz 
domain D — {x = {x,Xd) GM.'^ : Xd > jix)}- For such D the geometric notions 
of Theorem 2 become more explicit as we will see below. We note that D is 
((2V1 + A2)-i,r)-fat for ah r > [16, Remark 1]. For x = {x,Xd) G D and 



a/2 
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r > we define cc^''-' = (i, 7(i) + r). If a; is close to dD then x^^' dominates x in 
the direction of the last coordinate. We note that P^*"^' {jo > 1) ^ c > 0. Here 
c = c(c?, a, A). By Remark 3 and BHP, 

where C = C{a, d, A). By scaling, the Martin kernel with the pole at infinity for 
rD is a constant multiple of Mnix/r). By (44) we obtain 

P^rn >t) = P'-''^^{r,-u.^ > 1) ^ 1 A /I^.^L ' ^^^^ (^5) 

We note in passing that (45) agrees with (35) because r AIo {x^'^^ ) is increasing 
[15]. Or, in our previous notation we can take Ar{x, k, D) ~ 3;('"V(a;d-7(a;)) ^ 
substitute (45) into (1) so that for all < t < od and x,y £ D (in fact, by 
regularity, for a;, y G R'^) we have 

, c ( Md(x) \ , , f Mniv) 



Example 7. For circular cones V ([10]) we have 

Mv{x) = |x|^My(a;/|a;|) , x^Q, (46) 
where ^ /3 < a is a characteristics of the cone, see [1]. By [35, Lemma 3.3], 
My(a;)«^y(x)"/2|^|/3-'^/2^ X G M^ 

see also [10], [35]. Considering (44), by simple manipulations we obtain 

where C = C{\). By (1) and scahng we get (3). 

The interested reader may find more references on stable processes and Brow- 
nian motion in cones in [10]. Note that (46) holds for generalized open cones, 
i.e. open sets dW^ such that kV for all fc > [1]. 

Example 8. Let d — 1,2,... and V" = M"^' \ {a-d = 0}. This generalized cone 
is non-Lipschitz but it is (1/2, r)-fat for every r > 0. Let 1 < a < 2. From [1, 
Example 3.3] we have Mv{x) = \xd\°'~^ (the decay near a hyperplane is slower 
than near a half-space). We consider f = 1 in (36). We let Ai{x) = (i, Xd + 1/2) 
if > and Ai{x) = (i, a;^ — 1/2) otherwise. Thus, 



Mv(Ai(x)) (|a;rf| + l/2)^ 
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By (1) and scaling wc obtain the following analogue of (3), 



pit,x,y) V / V i^/" 

We note that V is the complement of a point if d = 1. 

If D is bounded and k > is fixed then D is not (k, r) at large scales 
r, and the asymptotics of the probability of survival is exponential. Indeed, 
for the fractional Laplacian with Dirichlet condition on D"^- we let Ai > be 
its first eigenvalue and 01 > the corresponding eigenfunction (normalized 
in L^{D,dx)) see [30]. The following approximation results from the intrinsic 
ultracontractivity of every bounded domain [30] : 

Here comparability constants depend on D and a (see also Proposition 1 below). 
Given that infinity is inaccessible [14] from bounded D, it is of considerable 
interest to understand the behavior of the heat kernel related to accessible and 
inaccessible points of D (see also [33] in this connection). 

In the remainder of the paper we will study C^'^ domains in more detail. We 
focus on unbounded domains, large times, and dependence of the comparability 
constants on global geometry of the domains. 

Example 1 and intrinsic ultracontractivity yield the following result. 

Lemma 6. There exist Ai ~ Xi{a,d) > and C = C{a,d) such that for all 
r > Q, t > and x € K'* we have 



a/2' 

I u T^l n ^ I \ 

1 A 



SB{Q.r){x) 

r A tV" 



Lemma 7. Let d > a, < r < R, W = B{0,r)U B'^{0, R). There is c = c{a,d) 
such that for all t > and x E M.'^ we have 

p-{rw>t)^c[^) 

Proof. By scaling we only need to consider r = 1 < R. By [5] we obtain 

^'<^«»''> - ~> - r((d-o)/2)rW2) i J^^"" » ' 

(compare (40)). Thus there is c = c{d, a) such that 



V rAfi/" J 



P^TBio.R)^^)>c>Q, \y\>2R. 
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Let a; e 5(0, 1). For t ^ 1 wc use (23) to obtain 



By (38), for t ^ 1 we even have 

P^(r.>t) ^ P^(r.(o.)>t)-lA(^f^2^) 



□ 



The C^'^ condition at a given scale fails to determine the fatness of D at larger 
scales and, consequently, the exact asymptotics of the survival probability. The 
following is a substitute. 



Proposition 1. If D is C^'^ at some scale r > then 



1 A 



5d{x) 
r A t^/° 



^ P^{td >t)^C 



1 A 



Soix) 



' A ti/c 



for all t > and x E R^. Here C = C{a,d) and Ai = Xi{a,d). 
If also d > a and diam(Z)'^) < oo, t/ien for all t > and x E R'^ , 



P^'iTD >t)^C- 



diam(_D'^ 



1 A 



Spjx) 
r A t^l°' 



a/2 



(49) 



(50) 



Proof. Consider x E D, <Z D and B{x" ,r) C D"^ of Remark 9. Clearly, 

tB:, ^ Tb(x",t), thus 

^"(tb. > i) P^'iro >t)^ P"(Tb(,",,) > i) . 
Lemma 6 yields the estimate 



1 A 



[ry 6d{x)) At^/"" 



a/2 



^ P%TD>t) 



and 



P%TD >t)i^C 



1 A 



Sd{x) 

r A tV" 



q/2' 



which simplifies to (49) a.s Sd {x) > r yields fo(a;)/[(r V (5d(x)) A i^/"] ^ 1. 
To prove (50), we consider p ~ diam(L''^) ^ 2r, the center, say xq, of i?^;, and 



K. Bogdan, T. Grzywny and M. Ryznar/Heat kernel estimates 



20 



W := BxU B''{xo, p + r\/ 5d{x)) C D. By Lemma 7 and Remark 9, 



^ C 



p + rV 5d{x) 



1 A 



Sd (x) 



[rV Sd{x)) At^^° 



a/2' 



1 A 



Spjx) 



a/2' 



□ 



In view of Theorem 1, (49) mildly strengthens [19, Theorem 1.1 (i)] (i.e. (2) 
above). We also get the following result. 

Theorem 3. Let d > a. If D is C^'^ at scale r and diam(L''^) < oo then 



c- 



diam(_D'^) 



2a 



< 



PD{t,x,y) 



1 A 



a/2' 



Pit,x,y) 



1 A 



a/2' 



for allt>0 and x,y e R"^. Here C = C{a,d). 
Proof. The result follows from (50) and Corollary 1. 

A similar result (with less control of the constants) is given in [22] ^ . 

Remark 10. We consider the recurrent case a^rf^l. IfDcRis the comple 
ment of a finite union of bounded closed intervals, then 

t>Q, x£W^ , if a > 1 , 
if a = 1 , 



□ 



P {TD>t) « 1 A ^i_i/^^^i/2 



s c- log(l + 5_d(x)^/2) , 



log(l+tl/2) 

where C = C{D,a). The estimates follow easily from Examples 2 and 3. 

Corollary 2. If D d W is the complement to a finite union of bounded closed 
intervals then C = C{D, a) exists such that for all t > and x, y G R, 



PDit,x,y) c 
P{t,x,y) 



1 A 



(5i3(.t)"-i A(5i5(x)"/2 



1 A 



Soiy)"-^ ASoiy) 



a/2 



for a > 1, while for a — I we have 



PDit,x,y) c 
P{t,x,y) 



1 A 



\og{l + 5D{xm 
log(l + ii/2) 



1 A 



\og{l + 5n{yy/'-) 
log(l + ti/2) 



-•^ paper [22] appeared on arXiv after the first draft [11] of the present paper 
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